Radial and angular rotons in trapped dipolar gases 
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We study Bose-Einstein condensates with purely dipolar interactions in oblate (pancake) traps. 
We find that the condensate always becomes unstable to collapse when the number of particles is 
sufficiently large. We analyze the instability, and find that it is the trapped-gas analogue of the 
"roton-maxon" instability previously reported for a gas that is unconfined in two dimensions. In 
addition, we find that under certain circumstances, the condensate wave function attains a biconcave 
shape, with its maximum density away from the center of the gas. These biconcave condensates 
become unstable due to azimuthl excitation - an angular roton. 

PACS numbers: 03.75Kk,03.75.Hh 
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The realization of a Bose-Einstein condensate (BEC) 
of 52 Cr Uli marks a major development in degen- 
erate quantum gases. The inter-particle interaction via 
magnetic dipoles in this BEC is much larger than that 
in alkali atoms, and leads to an observable change in 
the shape of the condensate. In these experiments, the 
dipolar interaction energy was about 15% of the short 
range interaction energy given by the scattering length. 
It may be expected that advances in utilizing Feshbach 
resonances |4| to tune the scattering length to zero, will 
lead in the near future to condensates dominated by pure 
dipolar interactions. A rich array of intriguing phenom- 
ena have already been pred icted to occur in such a gas 
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It is widely agreed that the structure and stability of 
a dipolar gas will depend critically on the interplay be- 
tween three factors: the anisotropic dipole-dipole interac- 
tion between particles; the anisotropic external potential 
used to confine the particles; and the number of parti- 
cles. Still, the details of this interplay remain thus far 
fragmentary. We take the standard scenario, where the 
dipoles are assumed to be polarized along a field whose 
direction defines the laboratory z axis. They are confined 
by a harmonic trap with cylindrical symmetry about this 
axis, defined by radial and axial angular frequencies u p 
and u z . It has been realized that, because the dipoles 
would prefer to align in a head-to-tail orientation, the 
gas elongates in the z direction and need not possess the 
same aspect ratio A = lo z /uj p as the trap HQ- If this 
is so, then ultimately the attractive interaction between 
dipoles should become so large as to initiate a macro- 
scopic collapse as the number of dipoles is increased. In 
fact, in the large-number limit, the gas approaches the 



Thomas-Fermi reg ime, in which a dipolar gas is inher- 
ently unstable |14| . 

The circumstances under which this collapse occurs has 
been the subject of much discussion. It was argued in 
Ref. that for oblate traps with A > 5.4, the gas is com- 
pletely stable, since the trap can overcome the tendency 
of the gas to elongate along z. However, Ref. [8j found 
unstable condensates even for A = 7. Taking it a step 
further, Ref. |l2| found instabilities even in the A — > oo 
limit, by treating a gas that is confined in the z direc- 
tion, but completely free to move in the transverse plane. 
The instability was refereed to as a "roton" instability, 
although the microscopic physics is quite different from 
the roton-maxon physics in superfluid He. Ref. E3 ex- 
amined the quasi-two-dimensional limit of this problem. 
Roton-like dispersion curves and supersolid-like phases 
had earlier been identified in quasi-lD gases, where the 
dipole dipole interaction is induced by strong laser light 

HEHillii. 

In any event, the details of the stability of a pure dipo- 
lar gas remain an open question. In this Letter, we re- 
visit this issue, applying a new efficient algorithm |l9| . 
We determine that at any finite value of A, the conden- 
sate will decay if sufficiently many dipoles are added. The 
route toward instability will be shown to exhibit physics 
similar to the roton mechanism, but adapted to a finite 
trap environment. Along the way, we demonstrate that 
stable condensates can have an unusual density profile, 
where the maximum of density is not in the center of the 
gas. 

The dynamics of the condensate wave function ip(r,t) 
is described by the time-dependent Gross-Pitaevskii 
(GP) equation: 
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FIG. 1: (Color online) Stability diagram of a dipolar conden- 
sate in a trap, as a function of the trap aspect ratio A = uj z /uo p 
and the dipolar interaction parameter D. Shadowed areas are 
stable, while white unstable against collapse. In the darker, 
isolated areas, we find biconcave condensates (illustrated with 
iso-density surface plot at the top right corner) whose maxi- 
mal density is not at the center. The contours in the biconcave 
regions indicate the ratio of the central density to the max- 
imal density, with darker areas having a smaller ratio. The 
contour intervals are 10%, and the minimum ratio obtained 
is 70%. 

where M is the particle mass and the wave function ip 
is normalized to unit norm. The coupling constant for 
the short-range interaction is proportional to the scat- 
tering length a, although we set a = in the follow- 
ing. The dipole-dipole interaction is given by Vd(r) = 
cP(l — 3cos 2 #)/r 3 , with being the angle between the 
vector r and the z axis. We also define a dimensionless 
dipolar interaction parameter, D = (N — l)Md 2 /(h 2 aho), 
where = J denotes the transverse harmonic os- 
cillator length. It is convenient to think of increasing D 
as equivalent to increasing the number of dipoles, for a 
fixed trap geometry. 

A stable condensate exists when there is a stationary 
solution of Eq. which is stable to small perturbations. 
Based on this criterion, we have constructed the stabil- 
ity diagram in Fig. ^ as a function of D and the trap 
aspect ratio A = lu z /lu p . In the figure, the shaded and 
white areas denote parameter ranges for which the con- 
densate is stable or unstable, respectively. In general, the 
more pancake-like the trap becomes (larger A), the more 
dipoles are required to make the condensate unstable. 

Remarkably, there appear regions in parameter space 



where the condensate obtains its maximum density away 
from the center of the trap. These are the darker shaded 
areas in Fig. ^ The local minimum of the density in 
the center gives the condensate a biconcave shape, re- 
sembling that of a red blood cell (a surface of constant 
density is illustrated at the top left corner of Fig. 
A density contour plot of such biconcave condensate is 
shown in Fig. |2{IIa), with the parameters A = 7 and 
D = 30.8. These structures appear in isolated regions of 
the parameter space, and in particular, only for certain 
aspect ratios in the vicinity of A w 7, 11, 15, 19.... There 
seems to be a repeated pattern which probably continues 
to larger values of A (although this was not calculated). 
However, the parameter space area of these regions be- 
comes increasingly small with larger A, and is already 
negligible for the fourth region (not shown in Fig. 
In between the biconcave regions, we find "normal" con- 
densates with maximum density in the center, Fig.[2IIa). 
Their density profiles in fact have a fairly sharp peak in 
the center, as compared to a Gaussian shape. 

To verify the existence of biconcave structures, we solve 
Eq. J3> numerically both with a 3D algorithm an d 
with a 2D algorithm that exploits the cylindrical symme- 
try [lij , carefully converging both the grid size and reso- 
lution. We have also carried out a variational calculation, 
in which the condensate wave function is taken to be a 
linear combination of two harmonic oscillator wave func- 
tions. The first is a simple Gaussian, and the second is 
the same Gaussian multiplied by (H2(x) + i?2(y)), where 
Hi is the Hermite polynomial of order 2, and (x, y) are 
the coordinates perpendicular to the trap axis. Minimiz- 
ing the variational energy analytically, we find biconcave 
solutions with a large component of the second wave func- 
tion, similar in shape to the numerical ones. The exact 
parameters for appearance of variational biconcave con- 
densates are somewhat different from the numerical ones, 
and the variational calculation gives a continuous param- 
eter region with biconcave condensates, rather than iso- 
lated regions as we find numerically. This is probably due 
to the oversimplified nature of the variational ansatz. 

It is perhaps not too surprising that a dipolar conden- 
sate could form with the density pushed to the outer rim. 
After all, the dipoles exert long-range, repulsive forces on 
one another, and flee to the surface, much as free charges 
in a conducting material do. What is somewhat more 
mysterious is why the dipoles do not always behave this 
way, but rather only for certain well-delineated islands in 
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FIG. 2: (Color online) (la) Contour plot of the ground state density for a dipolar condensate in a trap with aspect ratio A = 9 
and dipolar interaction strength D = 49.5 (close to instability).. (Ha) Condensate density contour plot for A = 7 and D = 30.8. 
Note the unusual shape, with the maximum of the density occurring away from the center of the trap, along a ring (the plot is a 
vertical cross-section through the cylindrically symmetric condensate). The 3-dimensional shape is biconcave (inset of Fig.0. 
(Ib) Excitation spectrum as a function of D for A = 9, and (lib) for \=7. Solid lines are excitation with angular quantum 
number m — 0, dotted lines m — 1, dashed lines m = 2 and dot-dashed lines m = 3. (Ic) The density perturbation (see text)of 
the lowest excitation mode for the ground state of (la), and the corresponding excitation (lie) for the ground states of (Ha). 
The heavy lines are nodal lines (actually surfaces in 3D) of the excitation. 



the stability diagram. 

Condensates with normal and biconcave shapes decay 
quite differently as the instability boundary is crossed. 
To examine the mechanism causing the collapse of con- 
densates, we have studied collective excitations in the 
Bogoliubov-De Gennes (BdG) approximation. Let us 
first consider the normal case, where the condensate den- 
sity is a maximum in the center of the trap. An example 
of this density profile is shown in Fig|2Ja. The excita- 
tion spectrum for this case is shown in Figl^Ib). The 
instability of the condensate at D — 50.0 coincides with 
an excited state mode "going soft," i.e., tending to zero 
excitation frequency. In this case the excited state has 
zero projection of angular momentum around the z axis, 
m = 0. The excitation consists of a radial nodal pat- 
tern, as we show in Fig.EJlc) by plotting contours of the 
eigenmode density Sn = u + v, where u, v are the usual 
BdG functions [19| . Thus when this condensate becomes 
unstable, the instability is due to a modulation of the 
condensate density in the radial direction. 

A similar modulation was already considered in the ro- 
ton instability studied for the infinite-pancake trap |l2j |. 
In that case, there is only a trap potential in the z di- 
rection, and motion is free in the (x, y) plane, so that 
the excitations have a definite transverse momentum q 
in this plane. The instability is expected to occur when 
this momentum approaches the value q « h/l z , where 
l z = ^Jh/Mijj z is the axial confinement length. In this 
case, excitations in the quasi-2D trap begin to acquire a 



3D character and the interparticle repulsion is reduced. 
For sufficiently many dipoles, the attraction can again 
overwhelm the condensate and initiate a collapse. Here 
the instability is a density-wave modulation sin(fcir) with 
characteristic wavelength 2-7r/fc ~ 2nl z . 

For a gas confined in all three dimensions, the same 
physics appears. For aspect ratio A = 9, the axial con- 
finement length is l z = a^o/3 (recall that the radial oscil- 
lator length defines our unit of length) . Thus the excited 
state should posses density- wave oscillations with a wave- 
length 2nl z ~ 2aho, or nodes at intervals ~ laho- The 
excitation plotted in Fig. &Ic) indeed has radial nodes 
on approximately this scale. In this sense, the confined 
gas exhibits what might be termed a "radial roton" in- 
stability. 

By contrast, for trap aspect ratios such as A = 7, where 
biconcave condensates appear [Fig.^IIa)], the lowest ex- 
citation mode near the instability boundary has angular 
momentum projection m > 0. For example, m — 3 for 
the case shown [Fig. Efllb)], meaning that the excited 
mode exhibits an azimuthal dependence proportional to 
sin(30). Therefore, with increasing number of particles, 
these condensates collapse due to density modulations in 
the angular coordinate, which spontaneously break the 
cylindrical symmetry Biconcave condensates thus de- 
cay to a kind of "angular roton" instability in the trap. 
[There is also one radial nodal surface, as seen in the 
(p,z) cross-section in Fig. [2{IIc).]. The appearance of 
an angular roton when the ground state has biconcave 
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(a) Zoom-in of Fig. 1 (b) Variational phase diagram 
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FIG. 3: (Color online) (a) Detail of Fig. showing double 
stability regions. The coarseness of the contour line is due 
to the resolution of our sampling of the parameters space, 
(b) Part of a stability diagram using a Gaussian variational 
method (there is no area corresponding to biconcave conden- 
sates, since these cannot be described by a Gaussian). 

shape may be understood in light of the fact that the 
maximum density of such a condensate lies a long a ring. 
The instability may thus also be described as buckling 
along this ring. 

In general, the number of the nodal surfaces of the 
soft mode near the instability increases with A: for radial 
rotons, from 3 for A = 9 to 7 for A = 17. For angular 
rotons, we find, for A = 6.6, m = 2 and one radial node; 
for A = 7, m = 3 and one radial node; and for A = 12, 
rn = 2 and 4 radial nodes. 

Finally, we point out an unusual feature in the stabil- 
ity diagram, shown in Fig. [21 Namely, for a small range 
of aspect ratios 6.36 < A < 6.5, the boundary between 
stable and unstable condensates bends back on itself, so 
that there are two distinct stability regions as D is var- 
ied. This feature is also reproduced qualitatively by a 
variational estimate using a single Gaussian |3fb)]. Ac- 
cording to the variational solution, we find that as D 
increases, the kinetic energy plays a role in de- stabilizing 
the condensate. Kinetic energy leads to a tendency to 
expand, but this tendency will be greater in the more 
tightly confined z direction. In this case, more dipoles 
can pile up along the axis, where the attraction between 
dipoles leads to instability. As D increases further, the 
dipole-dipole interaction becomes more important, and 
the effects of kinetic energy are mitigated. Finally, in 
the large-D limit, the Thomas-Fermi regime is reached, 
and complete instability occurs 0] . Evidently this dual- 
stability region depends an a delicate balance between 
the aspect ratio of the trap and the dipole-dipole interac- 
tion, since it occurs only over a small region of parameter 
space. 

In conclusion, we have described the stability dia- 
gram for dipolar condensates in pancake traps, found 
new structured, biconcave condensates, and an excitation 



spectrum which exhibits a discrete roton-maxon. The 
soft mode of biconcave condensates near collapse has an- 
gular momentum quantum number m > 0, i.e giving rise 
to azimuthal oscillations. This exotic behavior is another 
demonstration of the richness of the physics of dipolar 
condensates. 
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